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Recent progress in understanding general properties of high energy scattering near the unitarity limit, where high 
density gluon components of the wavefunction are dominant, is reviewed. The similarity of the QCD problem and 
that of reaction-diffusion processes in statistical physics is emphasized. The energy dependence of the saturation 
momentum and the status of geometric scaling are discussed. 

1 Introduction 

This lecture" deals with high density QCD where cross sections are near their unitarity limit (saturation). 
This is an area where there has been much new understanding and progress recently. Indeed, thre have been 
real breakthroughs in understanding between the time this lecture was originally given and the present time. 
The very recent developments will be alluded to, but not developed, somewhat later on. Conceptually, I 
believe we are near to achieving a rather complete understanding of how unitarity limits are reached in hard 
QCD processes at asymptotic energies. 

The lecture starts with a brief discussion of high energy soft scattering where a reasonable phenomenology 
exists but where there is not yet a good connection to QCD. Then the general issues of hard scattering at high 
energies are introduced with a focus on the idea of strong gluon fields and high gluon occupation numbers. 
The dipole picture of hard scattering is briefly reviewed and the Balitsky JIMWLK and Kovchegov equations 
obtained. Next a picture of the solution of the Kovchegov equation is given and the Y-dependence of the 
saturation momentum given. Then the problems with the mean field picture exemplified by the Kovchegov 
equation are given, and the essential non mean field aspects are outlined. This leads to an interesting rela- 
tionship between concepts entering high energy scattering in QCD and similar ideas currently under study 
concerning reaction-diffusion processes in statistical physics. 

2 High Energy "soft" Cross Sections 

Over a wide range of energies one can write, say for proton-proton collisions, ^'^ 



a*?*, ~ cis° °7 + C2/^A. (1) 

What is surprising about this result is that the form does not change significantly between fixed target energies 
and the highest energies at the Fermilab collider. One might have expected to see unitarity limitations more 
apparent in cr, however, one has to work harder to see that in fact unitarity limits are being reached even 
though the rate of growth of the total cross section has not changed. 

To see unitarity limits easily it is important to introduce an impact parameter dependence.^ Then 

J tot 

^ = 2(1-5(6,5)) (2) 



da"' 



= l-S' (3) 



= {l-SY (4) 
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where I have assumed that the S'-matrix is real for this simphfied discussion. Then the S-matrix can be 
obtained from proton-proton elastic scattering as 

with S schematically illustrated in Fig.l. At higher and higher energies S*^, the probability of not having 
an inelastic reaction, becomes very small at larger and larger impact parameters. When one integrates over 
impact parameters the form (1) mysteriously emerges which does not show the strong unitarity limits, and 
the energy independence of S for small b, exhibted in Fig.l. 




Figure 1. 



3 Hard Scattering at High Energies 

Consider the scattering of a small dipole, of transverse size x±, on a target at impact parameter b. The target 
can be a physical hadron or another dipole. When x± is very small'* 



2 Tra dxG{b,Q'^ =4:/xl) 

Sib, x^) = 1 - (6) 

where the energy (Bjorken-x) dependence has been suppressed. The .S-matrix only becomes significantly 
different from 1 when the target gluon distribution becomes of size 

dxG 



<PbdPk±_ a 
corresponding to an occupation number for gluons 



The scale at which gluon occupation becomes as large as in (8) and, hence, at which scattering is strong is 
called the saturation momentum.^ It could be more precisely defined by requiring, say. 



S{x^=2/Qs,b,Y) = l/2 (9) 

which equation gives Qs = Qs{b, Y) with Y the rapidity of the reaction. In terms of the strength of fields in 
the target (8) corresponds to ^4 l/gi at the scale of the dipole. 
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The key questions which arise and which wih be considered in much more detail in what follows are: (i) 
What is the F— dependence of Qs? (ii) What are the x± and F— dependences of T = 1 — 5? We recall that 
geometric scaling^ corresponds to 

T = T{xlQlb). (10) 

4 The Dipole Picture of Hard Scattering;^ the Balitsky-JIMWLK^-i° and Kovchegov^^ 
Equations 

Consider the scattering of a dipole on a target where the target carries almost all the momentum and is 
right-moving. The left-moving dipole has little momentum and servos as a probe of the highly evolved target. 
Now let F ^ y -f- dY and put the possible additional evolution into the probe which is a simple object. The 
initial dipole probe having a quark line at Xq and an antiquark line at x^ may have a virtual correction in 
the interval dY, in which case it remains an elementary dipole, or a gluon may now be present at X2 thus 
giving two dipoles. The two-dipole term consists of a dipole having a quark at and an antiquark at X2 and 
another dipole with a quark at and an antiquark at x^ . This gives the rapidity dependence of the scattering 
/S— matrix as 

dS{xoi,Y) _ aN, fd^ X2XQ1 
dY ~2^J 

[S^^\X02,X12,Y)-S{X0UY)] (11) 

where the S^'^^— term is the 5— matrix for scattering the two dipole term on the target and the 5— term is the 

virtual term. The factor ^Qid p jg ^j^^ element of probability for emitting a gluon, either real or virtual, 

^02^12 

and xoi = |^ — x^l, etc. Eq.ll is illustrated in Fig. 2. 




Figure 2. 



Eq.ll is part of the Balitsky-JIMWLK hierarchy.^~^° The difficulty with (11) is that one also needs an 

equation for S^"^"^ which will couple to a thrce-dipole scattering term, with the coupling to higher dipole terms 
never terminating. Kovchegov,^^ motivated by scattering off a large nucleus, suggested approximating 5*'^' by 

S^''\xo2,xi2,Y) = 5(xo2,F)5(xi2,r) (12) 

where all the 5— matrices are to be understood as referring to a given (but unspecified) impact parameter. 
Eq.l2 is a type of mean field approximation and leads to 

dS{xoi,Y) dNc f d'^X2xli 



• [S{x„2,Y)S{x,2,Y) - S{xoi,Y)]. (13) 

This equation, the Kovchegov equation, is now a precise equation for the scattering matrix and has a unique 
solution once the initial condition is fixed. Introducing T = 1 — S one may write (13) as 
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-^T{x,„Y) = ^J -^[T{xo2,Y) 

+T{xi2, Y) - T{xoi, Y) - T{xo2, Y)T{xi2, Y)] (14) 

with the first three terms on the right-hand side of (14) corresponding to the BFKL^^'^'' equation and the 
final term giving the nonhnear behavior which guarantees unitarity for the overall scattering. While (13) and 
(14) are not exact equations they are extremely useful as semi-realistic equations having BFKL behavior and 
obeying unitarity constraints. This is probably the best one can do in terms of a definite equation for a single 
function. The Balitsky-JIMWLK equations keep correlations better than the Kovchegov equation, but even 
they lack some of the essential ingredients necessary to impose unitarity in a genuinely realistic way.^^ 



5 The Solution (and Picture) of Evolution in the Kovchegov Equation 

While the Kovchegov equation is not an exact equation it does lead to interesting results some of which 
appear to be quite general and others which may be reasonably accurate in a limited energy region. One of 
the strongest results coming from the Kovchegov equation concerns the saturation momentun, Qg, introduced 
in (9). One finds5'i5-i8-i7 



2 _ 2aN, x(Ao) y 3£n{aY) 



1 - Ao 2(1 - Ao) 



1- Ao 



(15) 



inQl/A^ = J^^Y-cY^/^ (16) 
y TTO 1 — Ao 

Eq.(15) refers to a fixed coupling Kovchegov equation while (16) is for running coupling. The /x^ in (15) is 
some (scheme dependent) scale while the c in (16) is known in terms of the first zero of an Airy function. x('^) 
is the BFKL eigenvalue function with Aq satisfying 



x(Ao) = -(l-Ao)x'(Ao). (17) 

The first terms in (15) and (16), those involving the ratio ji^, were obtained long ago by Gribov, 
Levin and Ryskin^ from BFKL evolution. The full results were obtained in Ref.l7 using a cutoff to simulate 

the nonlinear terms in the Kovchegov equation. Munier and Peschanski^^ then obtained these full results 
using the general theory of nonlinear euqations in the FKPP^^'^° class. The FKPP equation is a "simple" 
one-dimensional equation 

dP(xx,t) d^P „ „o 

which has many of the general features contained in the more complicated Kovchegov equation. As we shall 
see a little later on, the leading asymptotic terms in (15) and (16) are correct although the exact results are 
different in the subleading asymptotic behaviors. 

In addition to the energy (rapidity) dependence of the saturation momentum the Kovchegov equation 
leads to a scaling property of the scattering amplitude which reads^ 

T{xouY) ^ ci[Q^(F)a;^i]i-^«[^n— V +C2] (19) 
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valid when Q^Xqi < 1 but where Q^x^i is not too small. The region where the Kovchegov equation gives 
scaling is quite a large region. The experimental data seem to support this scaling. Unfortunately, scaling is 
not correct asymptotically in a full treatment in QCD, but it may be true that there is an intermediate energy 
region where the Kovchegov equation is approximately valid giving (19) in the current HERA energy regime. 
This is an important point on which there is little understanding at the moment. 

In comparing results from the Kovchegov equation with experimental data one generally needs to use a 
model of the Golcc-Biernat Wiisthoff typc.^^'^''' How do the results hold up? Experimentally, -^InQ^/ ji^ ~ 
1/4. Eq.(15) gives a much larger value for -^inQl if one takes ~ 1/4, say. On the other hand, the 
running coupling answer given in (16) is much closer to experiment. Triantaiyllopoulos^^ has done a more 
complete analysis than that given in (16). He includes the higher order terms in the BFKL kernel as well as 
a resummation of coUinear effects along the lines of that developed by Ciafaloni, Colfcrai and Salam.^^ In the 
end he finds ^inQl ~ 1/4 - i in the HERA regime with the (pleasant) surprise being that higher order 
corrections and resummation do not change the result very much. 

6 Problems with the Mean Field Picture 

While the mean field picture contained in (13) and (14) can be expected to be accurate for scattering on a large 
nucleus over a limited energy regime these equations cannot in general be expected to be accurate. Numerical 
calculations by Salam^* sometime ago have shown that unitarity is reached in a very irregular way in high 
energy scattering. At one value of the impact parameter the scattering amlplitude, T, may be small while 
at another impact parameter T may be close to its unitarity limit. A similar phenomenon occurs at a fixed 
impact parameter when one considers an ensemble of possible wavefunctions. Some elements of the ensemble 
will have a a very small T while other elements of the ensemble may have T near one. Thus, an equation like 
(14), where such fluctuations are neglected, cannot be expected to lead to correct general results. What we 
shall see in the following sections is that the picture of fluctuations in the high energy hard scattering QCD 
problem (at a fixed impact parameter) closely resembles the fluctuations which are actively being studied by 
statistical physicists concerned with reaction-diffusion processes. 

7 Beyond the Kovchegov Equation; A First Try 

In trying to understand the domain of validity of the Kovchegov equation the authors of Ref.29 ran into the 
following dilemma. Imagine starting from an elementary dipolc of size xa s^t y = 0. Evolve the dipole, using 
BFKL evolution, to y = Y and ask what is the number density of dipoles of size xb &t Y. This is what one 
calls n{xA,XB,Y). n obeys a completeness relation 




(20) 



and since 



T{xa, Xb, F) ~ ca^x%n{xA, xb,Y) 
when x^ < x\, one can approximately rewrite (20) as 



(21) 




(22) 



where p = —£n{x^iJ,'^). 
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Suppose we imagine Xb being not too far from the saturation line, as shown in Fig. 3. Then in the 
approximation of using an absorptive cutoff to simulate the nonlinear term in the Kovchcgov equation the 
evolution described by Eqs.(20) and (22) is simply BFKL evolution with an absorptive boundary at pi{Y) as 
illustrated in Fig. 3. It is straightforward to see that the region of p which dominates the integral in (22) is 

p-Ps{Y/2)^VaV, (23) 
the normal diffusion region for BFKL evolution. But 



2 T{xA, X, Y/2) ~ e-(i-^°)(^-^^(^/2)) (24) 



1 

X 

in this region so that ^^(xa, x, Y/2) is typically extremely small in the integrand in (22). But if Xb is not far 
from the saturation line -^T{xa, xb- Y) is not small. This means that -^T{x^ XB^Y/2) must be typically very 
large. However, this latter fact is a violation of unitarity. Thus Ref. 29 concluded that while the Kovchegov 
equation gave unitarity limits on the overall evolution of a process it did not guarantee that intermediate 
stages of the evolution (going from y/2 to F in the above example) would obey unitarity limits. 

The simple procedure used in Ref.29 for "fixing" the problem was to limit the p— integration in (22) so that 
-^T(x, Xb, y/2) would never be larger than one. This was done by introducing a second boundary, P2(Y), to 
the right of which evolution was not allowed. The resulting problem was then that of simple BFKL evolution 
in the p — Y plane with the evolution being restricted to lie within the region pi < p < p2 as illustrated in 
Fig.3 

It is then not a difficult task to find the resulting saturation momentum to be^^ 

p^iY)=inQl/,^ = '-^f^Y 

TT i — Ao 

'''' 

where 



TT^x'XAo) 



2(Ap)2x(Ao) 



1 3 1 

Ap = —(nl/a^ H —Mn^ + const. (26) 

1 — Ao 1 — Aq 

Perhaps the bigger surprise is that geometric scaling is lost in the above approach with 



emerging rather than 



(27) 
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T = T{p-p,{Y)) (28) 

which follows from the Kovchegov equation. As we shall shortly see (25) is expected to be an exact asymptotic 
result while (27) represents a shortcoming of the "mean field" approach of Ref.29. 

8 Interpreting the Second Boundary; An Event by Event Picture of Evaluation 

Now, following Ref.30 and inspired by work on reaction-diffusion processes in statistical physics, we view 
evolution in a different way than we have presented it up to now. To have a specific example in mind consider 
the evolution of a "target" dipole which is elementary and of size xa sA, Y = 0. At rapidity Y we probe 
this evolved diple with an elementary dipole of size xb with the resulting scattering amplitude T{xa-,xb,Y) 
emerging. T{xa,xb,Y) is an average quantity because there are many possible states of the target, starting 
from an elementary dipole at y = 0, which contribute to T{xa, xb,Y). For example, in the procedure used by 
Salam^* to generate evolution of dipole wavefunctions using a Monte Carlo method one must sum over all the 
events generated at Y to obtain T{xa,xb-,Y). Now we wish to view evolution of a single event. For dipole 
evolution this evolution is governed by the same kernel which appears in (14). but now the kernel gives the 
probability that a dipole of size Xqi split into two dipoles of sizes Xq2 and 2:21. This is a stochastic process 
since what happens in a given "time" interval dy is not certain, but the various outcomes (splittings or no 
splitting) are given various probabilities, determined by the kernel of Eq.l4. Now denote by T{xA^XB^Y) the 
scattering amplitude for a dipole of size x b scattering on a particular component of the evolved state starting 
from xa at 2/ = 0. We call this a single event and note that 

<f{xA,XB,Y) >^T{XA,XB,Y). 

New view T{xa,Xb,Y) as a function oi p = —£n(p'^x'^). That is, for a particular evolution of the target 
starting from {xa) we make probes using projectile dipoles of various sizes, xb- The picture is as illustrated 
in Fig. 4. The picture is that of a travelling wave with velocity equal to ^^^|^- Since 

f{p, y) ~ co?n{xA ,XB,y), (29) 

and since n is discrete, the function T decreases rapidly to zero for p greater than p2 where T{p2,y) = c? ■ This 
intereprets the second boundary, discussed in the last section, as a cutoff reflecting the fact that for a given 
event the gluon occupation number cannot be less than one for any given "size" gluon. For Ps(y) < p < pi 
the evolution of T obeys the BFKL equation since gluon occupation is high and the stochastic part of the 
evolution is small. The approximation of Ref.29 is to neglect the stochastic part of the evolution completely 
and use BFKL evolution between the boundaries where T takes on values near one and near o? . 

What one gets out of viewing the process in an event by event way is the understanding that BFKL 
evolution is correct for the evolution of a single event away from T = 1 and T = o?. This gives one the 





1,2,3 are 
different events 



P 



Figure 5. 



understanding that the (25) should be an exact (asymptotic) result in QCD since the velocity for all travelling 
waves (events) is the same. 

What about the stochastic part of the evolution, the randomness inherent in the dipolc splitting when 
T ~ a^? The stochasticity is reflected in the fact that different events do not have exactly the same values of 
Ps, as illustrated in Fig. 5. The idea is that fluctuations in the tail of the distribution can lead, after their efi'ect 
diffuses into the main part of the travelling wave, to fluctuations in the saturation momentum from event to 
event. The strength of these fluctuations is given by 

=<pI> - < Ps >^ . (30) 
There have been elaborate numerical studies of in statistical models. The conclusion is that'^^'^^ 



(J = const. 



oY 



1/2 



(31) 



although, as yet, there is no analytical argument leading to this result. We presume that (31) also holds true 
for QCD evolution, at a definite impact parameter, and this naturally gives^° 



where the Ps{X) in (32) is the average, < PsiX) > of the saturation momentum coming from the various 
individual events. Eq.32 diff'ers from (27), where stochastic effects were ignored, and it also represents a 
violation of geometric scaling. 

The studies of statistical systems leading to (25) and (31) also suggest that these limits are reached 
rather slowly. '^"'^'■^^ Thus, while (25), (31) and (32) are likely exact results in QCD, although (25) needs some 
straightforward modifications to take running coupling effects into account, it is not yet clear that these effects 
can be seen at current energies. It is a very important task to better understand subasymptotic terms, both 
to connect with experiment and to see when asymptotia can be expected to arrive. Up to recently it has been 
thought that the Balitsky-JIMWLK equations are the correct equations for following amplitudes from weak 
to strong scattering in QCD as energy increases, and there has been an interesting numerical study^'^ of the 
Balitsky-JIMWLK equations. However, it now appears that the Balitsky-JIMWLK equations^'^*^ miss essential 
parts of the dynamics necessary to correctly impose unitarity on BFKL evolution. It will be interesting to 
see if a new equation emerges which corrects the deficiencies of Balitsky-JIMWLK and yet allow accurate 
numerical simulations so that the approach to asymptopia can be effectively studied. 
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